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>■ ' Abstract 

(^' 

In the framework of the Lindblad theory for open quantum systems, expres- 

^? , sions for the density operator, von Neumann entropy and effective temperature 

,__! I of the damped harmonic oscillator are obtained. The entropy for a state char- 

ON ' acterized by a Wigner distribution function which is Gaussian in form is found 

5^ . to depend only on the variance of the distribution function. We give a series of 

T^ \ inequalities, relating uncertainty to von Neumann entropy and linear entropy. 

II ■ We analyze the conditions for purity of states and show that for a special choice 

rH . of the diffusion coefficients, the correlated coherent states (squeezed coherent 

Cd I states) are the only states which remain pure all the time during the evolution 

of the considered system. These states are also the most stable under evolution 

in the presence of the environment and play an important role in the description 

of environment induced decoherence. 
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1 Introduction 

In the last years a large interest has arisen to study open systems, especially in quantum 
optics [H], quantum measurement theory |^ and in connection with decoherence and 
quantum to classical transition |^, ^, |], ||, 0|. The consistent description of open 
quantum systems was investigated by various authors ^ |^, |10|, 0, |12| (for a recent 



review see Ref. [|T3l). It is commonly understood [jTT], |lj] that dissipation in an 
open system results from microscopic reversible interactions between the observable 
system and the environment. Because dissipative processes imply irreversibility and, 
therefore, a preferred direction in time, it is generally thought that quantum dynamical 
semigroups are the basic tools to introduce dissipation in quantum mechanics. The 



most general form of the generators of such semigroups was given by Lindblad [15 
under the assumption that the evolution of the system is Markovian. This assumption 
is not generally valid, but it is considered a good approximation for various interesting 
physical models (for example, the quantum Brownian motion in a high temperature 
reservoir). In the Lindblad theory the density operator properties (Hermiticity, unit 
trace and positivity) are preserved by the master equation. 

Different master equations for the evolution of the density operator have been 
used to describe decoherence and quantum to classical transition ||^, ||, |^, |l^, |l^ |18|, 



T9| , P0| , |21| . In particular models, the density operator becomes diagonal in some 
basis, indicating that interference between the states in that basis is suppressed - 
this is the process of decoherence of density operators. In this case the dynamical 
variables corresponding to the diagonalizing basis are expected to evolve approximately 
classically. A different theoretical approach to the problem of quantum to classical 
transition is the decoherent histories approach [^, ^, ^, pl| . 

The Lindblad formalism has been studied for the case of the damped harmonic 
oscillator |jl3|, ^ |2^, ^ and applied to various physical phenomena, for instance, the 



damping of collective modes in deep inelastic collisions in nuclear physics |^. In |^| 
a family of master equations, constructed in the form of Lindblad generators, was pro- 
posed for local Ohmic quantum dissipation. In [^ the Lindblad master equation for 
the harmonic oscillator was transformed into Fokker-Planck equations for quasiproba- 
bility distributions and a comparative study was made for the Glauber P, antinormal 



ordering Q and Wigner W representations. In [^I], ^ the density matrix for the co- 
herent state representation and the Wigner distribution function subject to different 
types of initial conditions were obtained for the damped harmonic oscillator. A phase 
space representation for open quantum systems within the Lindblad theory was given 



m 



In the present paper we are also concerned with the observable system of a harmonic 
oscillator which interacts with an environment. In Sec. 2 we give the description of 
the Lindblad model for open quantum systems and write the master equation for the 
density operator of the damped harmonic oscillator. Generally the master equation 
gains considerably in clarity if it is represented in terms of the Wigner distribution 
function which satisfies the Fokker-Planck equation. In Sec. 3 we transform the 
master equation into the Fokker-Planck equation by means of the well-known methods 



111 , p^ , |35| , p6i |37| , p8| and analyze the evolution of the Wigner function of the density 
operator. The Wigner function is shown to have a Gaussian form. Then we derive 
a closed form of the density operator satisfying the master equation based on the 
Lindblad dynamics and describe the evolution of the considered system towards a final 
equilibrium state. By using the explicit form of the density operator, we calculate in 
Sec. 4 the von Neumann entropy and time dependent temperature of the quantum 
system in a state characterized by a Wigner distribution function which is Gaussian 
in form and analyze their temporal evolution towards the equilibrium values. Then 
we introduce the Wehrl entropy and present its basic properties, including the Wehrl- 
Lieb inequality and another inequality which gives an interesting relationship between 



the Wehrl entropy and von Neumann entropy. Following Halliwell and collaborators 
p9| , ^ , we give two other simple inequalites, relating uncertainty to von Neumann 
entropy, and von Neumann entropy to linear entropy. The dynamical behaviour of the 
Wehrl entropy is compared to the von Neumann entropy and Shannon information 
entropy (in some cases the two last entropies are equal to each other). The concept 
of the classical-like Wehrl entropy is a very informative measure describing the time 
evolution of a quantum system. The Wehrl entropy, first introduced as a classical 
entropy of a quantum state, can give additional insights into the dynamics of the 
system, as compared to other entropies. We also discuss about the introduction of 
entropy production for studying the stability of stationary states. In Sec. 5 we analyze 
under what conditions the open system can be described by a quantum pure state and 
show that for a special choice of the diffusion coefficients, the correlated coherent states 
(squeezed coherent states), taken as initial states, remain pure for all time during the 
evolution of the system. In some simple models of the damped harmonic oscillator 
in the framework of quantum statistical theory [^, ^, it was shown that the pure 
Glauber coherent states remain as those during the evolution and in all other cases, 
the oscillator immediately evolves into mixtures. In this respect we generalize this 
result, the results of other authors @], obtained by using different methods as well as 



our previous result of Ref. |^. In Sec. 6 we use the linear entropy for the description 



of the environment induced quantum decoherence phenomenon. We state that the 

correlated coherent states are the most stable under evolution in the presence of the 

environment and make the connection with the work done in this field by other authors 

I7| , pO| , p9| , ^ Q . Finally, a summary and concluding remarks are given in Sec. 7. 



2 Lindblad model for the damped quantum har- 
monic oscillator 



The simplest dynamics for an open system which describes an irreversible process is 
a semigroup of transformations introducing a preferred direction in time P, |l^, |15 



In Lindblad's axiomatic formalism of introducing dissipation in quantum mechanics, 
the usual von Neumann-Liouville equation ruling the time evolution of closed quan- 
tum systems is replaced by the following Markovian master equation for the density 



operator p{t) in the Schrodinger picture ||15|| : 

dMp) 



dt 



l{Mp))- 



Here, $t denotes the dynamical semigroup describing the irreversible time evolution of 
the open system in the Schrodinger representation and L the infinitesimal generator of 
the dynamical semigroup $t. The condition for the Lindblad theory is that the time 
scale considered for the open system (subsystem) should be very long compared to the 
relaxation time of the heat bath (external system), but shorter than the recurrence 
time of the total system assumed as a closed finite system [^, |4^. Using the structural 



theorem of Lindblad ||T5[ , which gives the most general form of the bounded, completely 
dissipative Liouville operator L, we obtain the explicit form of the most general time- 
homogeneous quantum mechanical Markovian master equation: 

^ = -UH,m] + ^uiy^m,v}] + fe-,p(t)v/]). (2) 



dt h' "^' '' 2h 



] 



Here H is the Hamiltonian of the system in the absence of environment. The operators 
Vj,Vj are bounded operators on the Hilbert space 7i of the Hamiltonian and they 
model the effect of the environment. 

We mention that the Markovian master equations found in the literature are of 
this form after some rearrangement of terms, even for unbounded Liouville operators. 



In this connection we assume that the general form of the master equation given by 
(0) is also valid for unbounded Liouville operators. 

We consider an open system consisting of a particle moving in a quadratic (har- 
monic oscillator) potential, coupled to an environment described by non-Hermitian 
Lindblad operators in (|^) which are a linear combination of position and momentum 
operators. We impose a simple condition on the operators H,Vj,Vj that they are 
functions of the basic observables q and p of the one-dimensional quantum mechanical 
system (with [g, p] = ihl, where / is the identity operator on 7i) of such kind that the 
obtained model is exactly solvable. A precise version of this last condition is that linear 
spaces spanned by first degree (respectively second degree) noncommutative polyno- 
mials in q and p are invariant to the action of the completely dissipative mapping L. 



This condition implies [25| that Vj are at most first degree polynomials in q and p and 
H is at most a second degree polynomial in q and p. Because in the linear space of 
the first degree polynomials in q and p the operators q and p give a basis, there exist 
only two C-linear independent operators Vi, V2 which can be written in the form 

Vj = ajp + bjq, j = 1,2, (3) 



with aj,bj complex numbers |2^. The constant term is omitted because its contribu- 
tion to the generator L is equivalent to terms in H linear in q and p which for simplicity 
are assumed to be zero. Then the harmonic oscillator Hamiltonian H is chosen of the 
general form 

H = Ho + ^{qp + pq), Hq = — / + ^-g'. (4) 



With these choices the Markovian master equation can be written [|T3], p^ : 



do % ^ % % 

-^ = -^[Ho,p] - ^(A + /i)[g,pp + pp] + ^(A-p)[p,pg + gp] 

-^[<i, [qJ]] - ^[P, \P,P]] + ^{[l \P,P]] + \P, [q,p]])- (5) 



Here we have used the notations: 

^ j=l,2 ^ i=l,2 

Dpg = Dyp = --Re ^ a*6j, A = -Im ^ a*6j, (6) 

i=l,2 i=l,2 

where Dgg, Dpp and Dpq are the diffusion coefficients and A is the friction constant. 



They satisfy the following fundamental constraints |13 



i)Dpp>0, ii)Z^,,>0, m)D,p,pDq,-Dl^>^. (7) 

The semigroup method is valid for the weak coupling regime, with the damping A 
typically obeying the inequality A <^ cuq, where uq is the lowest frequency typical of 
reversible motions. 

In the particular case when the asymptotic state is a Gibbs state 



Pg{oo) = e kT /"Tie ^^^ (8) 



these coefficients reduce to 



A + /i, ^ hu \ — u h ^ hu 

nmuj coth — — -, Dan = coth — — -, 

2 2kT' ^^ 2 mu 2kT' 



Drm = ^mtucoth— — -, Daa = coth — — -, Dr,„ = 0, (9) 



where T is the temperature of the thermal bath and the fundamental constraints are 
satisfied only if A > |yu|. 

The necessary and sufficient condition for L to be translationally invariant is A = /i 



||13| , ^ , |26| . Translation invariance means that [p,L{p)] = L{[p,p]). In the following 
general values for A and /i will be considered. 

In the hterature, master equations of the type (|) are encountered in concrete 
theoretical models for the description of different physical phenomena in quantum 
optics H, ^ ^, ^ , in treatments of the damping of collective modes in deep inelastic 



collisions of heavy ions [Q or in the quantum description of the dissipation for the 



one-dimensional harmonic oscillator |]TT], |I2|, ^, ^ . A classification of these equations 



whether they satisfy or not the fundamental constraints (|^, was given in ^7j. Density 
operators satisfying the Lindblad equation give statistical predictions in full agreement 
with experiment in a wide variety of situations. Lindblad type equations are also 
frequently used in studies of quantum decoherence |Q, |^, |T^, |2T|]. For example, in the 
much-studied quantum Brownian model |51, 52 1, the master equation is the Lindblad 



master equation with a single Lindblad operator 

V = {2D)-'/\q + 2^^Dp), (10) 

with D = h'^ /SrwykT (where 7 is the dissipation and T is the temperature of the 
environment) and H = Hg + (7/2){g, p}, where Hg is the distinguished subsystem 
Hamiltonian in the absence of the environment. The Lindblad master equation does 
not, in fact, completely agree with the master equations given in a number of papers on 
quantum Brownian motion. In particular, the master equation given by Caldeira and 
Leggett [^ does not contain the term [p, [p, p]] and is known to violate the positivity 
of the density operator on short time scales |S^. This difference is principial, but 
practically Lindblad and Caldeira-Leggett equations are similar for high temperature, 
when the Markovian approximation is only valid and the extra term is negligible, since 
its coefficient is proportional to T^^ |^ . 



In the following we denote by gaa the dispersion of the operator A, i.e. aAA =< 
A^ > - < A >'^, where < A >= aA = Tr(pi) and Trp = 1. By aAB = 1/2 < 
AB + BA > — < A >< B > we denote the correlation of operators A and B. 

From the master equation (M) we obtain the following equations of motion for the 



expectation values and variances of the coordinate and momentum [13, It 



^ = -(A-p)a,(t) + la,(t), (11) 

at rn 

dop{t) _ 



dt 



muo'a^{t)-{\ + ^i)ap{t) (12) 



and 



^^^ = -2(A - /i)(T,,(t) + ^o-p,(t) + 2D,„ (13) 

^ = -2(A + /x)(Tpp(t) - 2m^Vp,(t) + 2Dpp, (14) 

^^^ = -muj^aqqit) + -app(t) - 2A(Tp,(t) + 2/}^,. (15) 

LLL lit 

In the under damped case {u > ji) considered in this paper, with the notation Vt^ = 
co^-fi^, we obtain 0,^: 

aJt) = e-^\(cos nt + ^ sin nt) a J 0) + ^ sin nta JO)), (16) 

il mil 

2 
TflUJ Li 

ap{t) = e"^*( —sinntag{0) + (cos fit- ^ sinfit)ap(O)) (17) 

and crg(oo) = crp(oo) = 0. It is convenient to consider the vectors 



^ muaqqit) 
X{t) = app{t)/muj 

\ ^pqif) 

and 

2mujD 



D = I 2Dpp/muj I . (19) 

2Dr 



'pq 



With these notations the solutions for the variances can be written in the form [|T3], pB| : 

X{t) = (Te-^*T)(X(0) - X(oo)) + X(oo), (20) 

where the matrices T and K are given by 

-, / fi + iQ fi — iQ 2uj 

^-iVl ^ + iVl 2uj I , (21) 



2ifi , o 

—uj —u —2jj, 



2{X-iQ) 

K= \ 2{X + tn) I (22) 

2A 



and 

X(oo) = {TK~^T)D. (23) 

The formula (^) is remarkable because it gives a very simple connection between 
the asymptotic values (t -^ oo) of cXqqlt) , (Tpp{t) , apq{t) and the diffusion coefficients 

^qqi ^ppi pg ■ 

(m2(2A(A + /i) + uj^)Dqq + Dpp + 2m(A + fi)Dpq), (24) 

"t" LU" — fX" ) 
1 



{{mujfuj^Dqq + (2A(A - /i) + cj2)Dpp - 2mcj2(A - /i)Dp,), (25) 



-PPV--/ 2A(A2 + cu2 - /i2) 

o^P.(oo) = 2^A(A2 + a;2-/x2) ^~^^ + /i)(m^)2D,, + (A - ^i)Dpp + 2m(A2 - ^?)Dpq).{2Q) 

These relations show that the asymptotic values crqg(oo), crpp(oo), o"pq(oo) do not depend 
on the initial values crgg(O), crpp(O), crpq(O). In the considered underdamped case we have 

„-2At / &11 ^12 &13 



rp^-Ktrp 



&21 &22 ^23 I , (27) 

\ &31 ^32 &33 



2fi2 

where btj, i,j = 1, 2, 3 are time-dependent oscillating functions given by (3.78) in [p6| . 



3 Wigner distribution function and density opera- 
tor 

One useful way to study the consequences of the master equation (|^) for the density 
operator of the one-dimensional damped harmonic oscillator is to transform it into 
more familiar forms, such as the equations for the c-number quasiprobability distribu- 
tions Glauber P, antinormal ordering Q and Wigner W associated with the density 



operator |3y] . In this case the resulting differential equations of the Fokker-Planck type 



for the distribution functions can be solved by standard methods p4|, ^, ^, ^ ^ 



10 



employed in quantum optics and observables directly calculated as correlations of these 
distribution functions. 

In [^, ^, Q we have transformed the master equation (|^) for the density operator 



into the following Fokker-Planck equation satisfied by the Wigner distribution function 

W{q,p,t): 

dW V dW o dW d d 



dt m dq dp dq dp 

d'^W aw d'^W 

Qq2 PP Qp2 ^'^dpdq 



+^^^1^ + ^PPl^ + ^P.p-^- (28) 



Since the drift coefficients are linear in the variables q and p and the diffusion coeffi- 
cients are constant with respect to q and p, Eq. (|28|) describes an Ornstein-Uhlenbeck 



process [^. Following the method developed by Wang and Uhlenbeck [^, we solved 



in ||3T| this Fokker-Planck equation, subject to either the wave packet type or the 5- 
function type of initial conditions in the underdamped case {u > ^) of the harmonic 
oscillator. One gets two-dimensional Gaussian distributions with different variances. 
Wigner function allows us to compute the moments of q and p at any time t in terms 
of the initial moments. By computing the long time limits of these moments, the form 
of the long time limit of the Wigner function may be obtained, since it is completely 
determined by its moments. 

For an initial Gaussian Wigner function the solution of Eq. (pSD is 

1 



W{q,p,t) 



2n^ 



X exp{-^[appit)iq - a.it)^ + a,,{t){p - a^{t)f - 2a^,{t){q - a,(t)){p - ap(t))]},(29) 
where a{t) is the determinant of the dispersion (correlation) matrix M{t), 

a{t) = det M{t) = a,,{t)app{t) - aj^(t) (30) 

and 

M(t) = f^'''^5J| ^p.wy (31) 

11 



We see that the initial Wigner function remains Gaussian and therefore the property 
of positivity is preserved in time. When time t — ^ oo, (Tq(t) and o"p(t) vanish, all 
dependence on crg(O) and o"p(0) drops out of the exponentials in W and we obtain the 
steady state solution: 

Woo{q,p) = ^==exp{-— - — -[o-pp(oo)g^ + (Tgg(oo)/ - 2apq{oo)qp]}, (32) 

27r./(T oo 2a(oo) 



where a{oo) = 0-^^(00) a pp{oo) — 0-^^(00) and crgg(oo), crpp(oo), crpg(oo) are given by Eqs. 
(H) ^ (iD- All stationary solutions to the evolution equations obtained in the long 
time limit are possible as a result of a balance between the wave packet spreading 
induced by the Hamiltonian and the localizing effect of the Lindblad operators. 

We obtain now the explicit form of the density operator of the damped harmonic 
oscillator in the Lindblad theory by using a technique analogous to those applied in the 
description of quantum relaxation |^ ^ ^ ^. Namely, we apply, like in [^, ^ , 
the relation p = 27ihN {Ws{q, p)} , where Wg is the Wigner distribution function in the 



form of standard rule of association and N is the normal ordering operator ||3^, ^ 
which acting on the function Ws{q,p) moves p to the right of q. By the standard rule 
of association we mean the correspondence p'^q" -^ qnpm Iqq^^t/^qq^ functions of two 
classical variables {q,p) and functions of two quantum canonical operators (g, p). The 



Wigner distribution function (^), which is in the form of the Weyl rule of association 
||60| , can be transformed into the form of standard rule of association ||61[ by performing 
the operation: 

Ws{q,p) = exp{^th-^)W{q,p). (33) 

The normal ordering operator A^ can be applied upon the Wigner function Ws in 



Gaussian form by using McCoy's theorem ||35|, |59|. Following Jang [56], we obtained 



in [^ the following expression of the density operator: 



p[tj = —j= expj- In — — -2 cosh (1 + 



12 



X [a^pim - a,(t))2 + a,,it)ip - a,it)f - 2a,,it)[iq - a,(t))(p - a^it)) - 2^]]}, (34) 
where 

b = CXgq{t)app{t) - {cxpq{t) - l-f . (35) 

The density operator (Q) has a Gaussian form, as expected from the initial form of 
the Wigner distribution function. While the Wigner distribution is expressed in terms 
of real variables q and p, the density operator is a function of operators q and p. When 
time t -^ cx), the density operator tends to 

p(oo) = 



Aaioo) - h^ 



1 2Ja{oo) + h 

xexp{ . In ^ [app{oo)q + agg{oo)p - crpg{oo){qp + pq)]}. (36) 

2hJa{oo) 2Ja{oo) — % 



In the particular case (^, 

(Too ( oo ) = coth— — -, a„Joo) = coth— — -, cr„o(oo) = (37) 

*^^ ^ 2mcu 2kT' ^^^ ^ 2 2A;r' ^^^ ' ^ ' 

and the asymptotic state is a Gibbs state (1): 

Pc(oo) = 2=mh — exp{--(-/ + —,')}. (38) 

The Wigner function can be expressed as the Fourier transform of the off-diagonal 
matrix elements of the density operator in the coordinate representation ||62|| : 



W{q,p) = ^Jdy <q-y\p\q + y> e''^y/\ (39) 

Then < x\p\y > can be obtained from the inverse Fourier transform of the Wigner 
function: 

/t X ~\~ V 

dpexp{-p{x~y))W{p,^—), (40) 
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namely 

-^KpW - ^)(^ - »)' + 1^(^ - "Mi^ -y) + j'y.m;': - v)]m 

In the long time limit crg(t) = 0, crp(t) = and we have 

< Ap{^)\y >= (t^ — r^)^ 

2TXaqq[00) 

X exp[-^-4-T(x + y? - ^K.(oo) - ^4^)(^ - y? + ^^(^' - y%m 
4 Entropy and uncertainty 

The most natural measure of uncertainty of the quantum mechanical state is the 
entropy. Physically, entropy can be interpreted as a measure of the lack of knowl- 
edge (disorder) of the system. The von Neumann entropy measures deviations from 
pure state behaviour. For an isolated system the entropy is time independent due 
to the unitarity of the evolution operator. For open systems such as the damped 
harmonic oscillator, the evolution is not unitary and the entropy becomes time depen- 
dent. Denoting by p{t) the density operator of the damped harmonic oscillator in the 
Schrodinger picture, the von Neumann entropy S{t) is given by the expectation value 
of the logarithmic operator Inp (we put Boltzmann's constant k = 1) : 

S{t) = - < Inp(t) >= -Tr(p(t) Inp(t)). (43) 

Accordingly, the calculation of the entropy reduces to the problem of finding the ex- 
plicit form of the density operator. Using Eq. (34) we obtained in [3^ the following 
expression of the von Neumann entropy: 



^(t) = (z/ + l)ln(z/ + l)-z/lnz/, (44) 



where we denote hu = Jo-(t) — h/2. It is worth noting that the entropy depends 
only upon the variance of the Wigner distribution. When time t ^ oo, u tends 

14 



to s = Ja{oo)/h — 1/2 and the entropy relaxes to its equilibrium value S{oo) = 
{s + 1) ln(s + 1) — s In s. The expression (H3) is analogous to those previously obtained 



p3[| in the theory of quantum oscillator relaxation and for the description of a 



system of collective RPA phonons [|56|. It should also be noted that the expression 



has the same form as the entropy of a system of harmonic oscillators in thermal 



equilibrium. In the later case u represents the average of the number operator ^5 
Although the expression (|4^) for the entropy has a well-known form, the function 
z/ induces a specific behaviour of the entropy. From the expression of the variances 
(pop, (pTl), which appear in (|30D , we can see that the time dependence of the entropy 
is given by the damping factors exp(— 4At), exp(— 2At) and the oscillating function 
exp 2iQt. The complex oscillating factor exp 2iQt reduces to a function of the frequency 
u, namely exp2iujt, for /i ^ or if fx/Q ^ 1 (i.e. the frequency u is very large as 
compared to fi). 

The von Neumann entropy gives zero for all pure states p = p^, so it measures the 
purity of quantum states, being different from zero (in fact positive) only for mixed 
states. This entropy does not differentiate between various pure states. We have pure 
states for a = h^ /A, which is just the case when we have equality in the generalized 
uncertainty relation a > fi /A |^ (see Sec. 5). By calculating 



- = lu(l + -). (45) 

we see that S is increasing when u increases. But in general the entropy is not a 
monotonic function of time, because we cannot decide a priori about the sign of dS/dt. 
In fact, in the underdamped case considered in this paper, the entropy oscillates and 
only when t — > cxd it relaxes to its equilibrium value, which is independent of the initial 
state of the system. If we require that dS{t)/dt > 0, then the environment operators 
Vj must be Hermitian. This assertion was proved in Ref. [^] by a direct calculation 
and the use of the inequality (x — y) \og{x/y) > for positive x, y. 

15 



In the case of a thermal bath (|]), a time dependent effective temperature Te can 
be defined |56| , p7| , by noticing that when t -^ oo, u tends, according to (2.24) in Ref. 



?T| , to the average thermal phonon number < n >= {exp{huj / kT) — 1)^^. Thus z/ can 
be considered as giving the time evolution of the thermal phonon number, so that we 
can put in this case 

(exp — -l)-^ = z/. (46) 

From (146|) the effective temperature Tg can be expressed as 

K[ln(z/ + 1) — mz/J 
Accordingly, we can say that at time t the system is in thermal equilibrium at temper- 
ature Te. In terms of the effective temperature, the von Neumann entropy takes the 



form p2| 



■5W = ih * '"11 - '=M~)l (48) 

r.(exp|i-l) *^' 



- e 



As t increases, the effective temperature approaches thermal equilibrium with the bath. 

The von Neumann entropy of the density operator is often connected with uncer- 
tainty, decoherence and correlations of the distinguished system with its environment 
0, ^, |l^ |18|, ^. Zurek, Paz and Habib, for example, looked for classes of initial states 



which generate the least amount of entropy at time t. They regarded such states as 
the most stable under evolution in the presence of an environment. They argued that 
these states are coherent states, at least approximately. The coherent states are known 
to be as close as possible to classical states: quantum fluctuations are minimal in these 
states and equal to those of the vacuum. Anastopoulos and Halliwell have shown that 
it is really the correlated coherent states, rather than the ordinary coherent states 
which are the most stable |^. One of the reasons for looking at the von Neumann 

16 



entropy is that it is constant for unitary evolution, thus for open systems such as that 
considered here, it is mainly a measure of environmentally induced effects (see Sec. 6). 
Entropic uncertainty relations have been obtained by Bialynicki-Birula in Refs. 
pB|, 1^. These relations express restrictions imposed by quantum theory on probability 
distributions of canonically conjugate variables in terms of corresponding information 



entropies [3C]. Entropic uncertainty relations derived from phase space quasiprob- 
ability distributions and connection to the decay of quantum coherences have been 



discussed in Refs. [£8|. Further discussion on entropic uncertainty relations associated 



with the phase-number uncertainties can be found in Ref. |69 



Following 1^9], ^, we shall give the connection between uncertainty and the von 
Neumann entropy, by considering the phase space quasiprobability Qh distribution 



QHip,q) =< Oi\p\a >, 



(49) 



where 



< x\a >=< x\p, q >- 



1 



(2vr.,,)V^ 



exp( 



(x — q)" 



4s 



<?? 



+ -px) 



(50) 



are the standard coherent states, Sqq and Spp are the variances of the coherent states 
and SqqSpp = h'^ /A. It is normalized according to 

dpdq 



/^<2«<-'-^^ 



(51) 



Between the two quasiprobability distributions Qh and Wigner W there is the follow- 
ing relation [Q: 



Qnip^q) = 2 / dp'dq' exp 



(p - p'Y (q - q') 



n2' 



2s 2s 



W{q',p'). 



(52) 



Eq. (|52D is the Husimi distribution |7^. In contrast to the other two quasiprobability 
distributions Glauber P and Wigner W, the Qh representation is always a positive 
and well-behaved function [^, |71|. The distribution QniPiq) is therefore a Wigner 
function, Gaussian smeared over an ^-sized region of phase space. 
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An information-theoretic measure of the uncertainty in phase space contained in 
the distribution (0), is given by the Shannon information of the Qh distribution: 

I = -J^QHip,q)\nQH{p,q). (53) 



This is the so-called Wehrl entropy |7^ , defined in full analogy to the classical entropy 
in phase space. Wehrl considered this entropy as a classical analogue of the von 
Neumann entropy and it was introduced as a classical entropy of a quantum state. 
It gives additional insight into the system dynamics, as compared to other entropies. 
/ is large for spread out distributions and small for very concentrated ones. The 
uncertainty principle implies that a genuine phase space probability distribution in 
quantum mechanics cannot be arbitrarily peaked about a point in phase space. The 
information (^) should possess a lower bound and since the coherent states are the 
states most concentrated in phase space, we expect the lower bound to be the value of 
J on a coherent state [^ . The most important property of this classical-like entropy 
is the following inequality, which gives a lower bound on the Shannon information: 

/ > 1, (54) 

with equality if and only if the considered state p is a coherent state [^ . The Wehrl 
entropy is a good measure of the strength of the coherent component and it clearly dis- 
tinguishes coherent states, in other words it measures how close a given state is to the 
coherent states or how much coherence a given state has. Wehrl entropy is very sensi- 
tive to the phase space djTiamics (such as, e. g., spreading) of the Qh representation. 
It extracts from the Qh function essential information about the investigated system. 
The Wehrl entropy cannot be negative. This follows from the fact that < Qh{oi) < 1 
and from the normalization condition (|5T|). The Qh function can never be so concen- 
trated as to make / negative. On the contrary, classical distributions can be arbitrarily 
concentrated in phase space and classical entropies can take on negative values. They 
may even tend to — oo if the distributions tend to (5-functions. 



The Shannon information can be used ^^ to find a relation between the von 
Neumann entropy and the generahzed uncertainty measure 



a = aqqUpp - CTpg. 



(55) 



The Shannon information satisfies the inequahty 



J < ln(^(det(T(^))^ 



(56) 



where o'^^^ is the 2x2 covariance matrix of the distribution Qh{p, q) 0, 0- Equahty 
holds if and only if Qh{p, q) is a Gaussian. From (|52D or using the results of Ref. [^ 
one has 



det or('5») = (agq + Sqg) {app + Spp) - a' 



2 
pq- 



(57) 



There is an even stronger relationship 1)72 



/ > ^, 



(5^ 



which establishes a connection between the Wehrl entropy and the von Neumann 
entropy of a given state. Using the relations ( ^ 



(|58D , one obtains 



{cTqq + Sqq) {(Tpp + Spp) - alg > k^e^'^^ ^\ 



(59) 



Since Sgq is arbitrary and SggSpp = Tt? 1^-, we may minimize the left-hand side over it and 
we obtain the following connection between the uncertainty and entropy for a general 
mixed state p [^ : 



(v^ 



+ ?)'-ai > ^V(^-^). 



O'qqCTpp , . ^ pg 



(60) 



An analogous relation was obtained by Dodonov and Man'ko in |75 |. 

In the regime where quantum fiuctuations are more significant than thermal ones, 
it is appropriate to use the lower bound (|5^) rather than (|58| ) (since S goes to zero 
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if the state is pure) and this is formally achieved by setting S" = 1 in (|60D. One then 
deduces the usual uncertainty principle from ( |60D ||40[] . 

In the regime where thermal (or environmentally induced) fluctuations are domi- 
nant, one would expect aggapp ^ h^ /A and 5^1 and (^) then gives 

V^^^ > rl > e^. (61) 

Here A = y/a is the Wigner function area - a measure of the phase space area in which 
the Gaussian density matrix is localized and y/agqapp/Ti is the number of phase space 
cells occupied by the state [^ . The von Neumann entropy of a Gaussian is given by 



and for large A^ (|4^) gives 



S^\nj (62) 



and hence we have equality in (|61|) 



We can conclude that J is a useful measure of both quantum and thermal fluctua- 
tions. It possesses a lower bound expressing the effect of quantum fluctuations and is 
closely connected to entropy, which in turn is a measure of thermal fluctuations ||39|| . 



In some models the linear entropy is introduced as a measure of purity of states: 

Si = Tr(/3 - p2) = 1 - Trpl (63) 

Since Trp^ < 1, the linear entropy is positive and it becomes zero if the state is pure. 



When the linear entropy is increasing, the degree of purity is decreasing. In [|^ the 
following relation is deduced between the linear and von Neumann entropies: 

^/ = 1 - Trp2 < 1 - e--^. (64) 

Equality in Eq. (|6^) is reached for pure states, when S = Si = Q and for very mixed 
states, when S is very large and S*/ ~ 1. 

For analyzing the stability of stationary states, it could also be interesting to intro- 
duce the entropy production crs^p) |7^, defined for dynamical semigroups as follows 
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asip) = -f^S{Mp)\p%=o, peLliH), (65) 

where 

S{Mp)\p') = Tr($*(p) lnp°) + S($,(p)), (66) 

p^ G L^iTi) is a $t-invariant state, i.e. ^tP^ = P^, for any t > 0, $t(p) = p(t) and 
L\_{7i) is the Banach space of trace class operators on 7i. 

5 Purity of states 

Schrodinger ||78[ and Robertson [^ proved for any Hermitian operators A and B and 
for pure quantum states the foUowing generahzed uncertainty relation : 

cTAAcrBB -o-Ij,>-\< [A, B] > |2. (67) 

For the particular case of the operators of the coordinate q and momentum p the 
uncertainty relation (^) takes the form 

CTgqapp - al > — . (68) 
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This result was generalized for arbitrary operators (in general non-Hermitian) and 
for the most general case of mixed states in [Q. The inequality ( |68|) can also be 
represented in the following form: 

^^g^pp - 4(1 -r2) ' ^^^' 

where 

(70) 



v^ 



O'qqO'pp 



is the correlation coefficient. The equality in the relation (|68D is realized for a special 
class of pure states, called correlated coherent states [|^ or squeezed coherent states, 
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which are represented by Gaussian wave packets in the coordinate representation. 
These minimizing states, which generahze the Glauber coherent states, are eigenstates 
of an operator of the form |^^ : 



1 IT T] 

with real parameters r and r], \r\ < 1,1] = Jo\[q. Their normalized eigenf unctions, the 



correlated coherent states, have the form |64 

with a. a complex number. If we set r = and r] = (h/2mujY^'^, where m and u are the 
mass and respectively the frequency of the harmonic oscilator, the states (|72D become 
the usual Glauber coherent states. In Wigner representation, the states ( [72D have the 



form []64| : 



This is the classical normal distribution to give dispersion 

2 fi hr 

<ygq = V , (^pp = 4^2(1 _ ^2)' ^W = 2(l-r2)l/2 ^'^^^ 

and the correlation coefficient r. The Gaussian distribution ( [75D is the only positive 



Wigner distribution for a pure state ||8^. All other Wigner functions that describe 
pure states necessarily take on negative values for some values of q,p. 

In the case of the relation ( |67D the equality is generally obtained only for pure 
states [0. For any density matrix in the coordinate representation (normalized to 
unity) the following relation must be fulfilled: 

7 = Trp2 < 1. (75) 

The quantity 7 characterizes the degree of purity of the state. For pure states 7 = 1, 
for highly mixed states 7 <^ 1 and for weekly mixed states 1 — 7 ^ 1. In the language 
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of the Wigner function the condition (^) has the form: 

7 = 27r^ / W^{q, p)dqdp < 1. (76) 

Let us consider the most general mixed squeezed states described by the Wigner 
function of the generic Gaussian form with five real parameters: 

W{q,p) = - — ^exp{- — [app{q - aqf + aqq{p - Opf - 2apq{q - aq){p - dp)]}, (77) 
I'RsjG la 



where a is given by Eq. (pSf ). The Gaussian Wigner functions of this form correspond 
to the so-called quasi-free states on the C*-algebra of the canonical commutation re- 
lations, which is the most natural framework for a unified treatment of quantum and 
thermal fluctuations [Q. For Gaussian states of the form ( [77D the coefficient of purity 



7 is given by 



^ - h- <^'' 



Therefore, we reobtain from ( ffSD and ([78|) that a has to satisfy the Schrodinger- 
Robertson uncertainty relation ( |68D 

a>-. (79) 

This inequality must be fulfilled actually for any states, not only Gaussian. Any 
Gaussian pure state minimizes the relation (^). For a > fi /A the function ( [77[) 
corresponds to mixed quantum states, while in the case of the equality a = h^ /A it 
takes the form ( [73D corresponding to pure correlated coherent states (squeezed coherent 
states). 

We have seen in the preceding section that the degree of purity of a state can also 
be characterized by other quantities besides 7. The most usual one is the quantum 
(von Neumann) entropy. For quantum pure states this entropy is identically equal to 
zero. We remind [0, |5^ that for Gaussian states with the Wigner functions (^) the 
entropy can be expressed through a only, in the form (0). 
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In the Lindblad model for the damped harmonic oscillator, the relation (|^ is 
a necessary condition for the generalized uncertainty inequality ( ]68| ) to be fulfilled 
| |T^ , [2^ , P^ . By using the fact that the linear positive mapping defined by A — * Tr(/3A) 
is completely positive, the following inequahty was obtained in Refs. |]T^, PB[ : 

Dpp(Jqq{t) + Dqqapp{t) - 2Dpq(rpg{t) > — . (80) 

We have found in [^| that this inequality, which must be valid for all values of t, is 



equivalent with the generalized uncertainty inequality (|68D at any time t 



a,,(t)app(t)-<(t)>-, (81) 

if the initial values crgg(O), o"pp(0) and <Jpq{0) for t = satisfy this inequality. If the 
initial state is the ground state of the harmonic oscillator, then 

h , . mhuo , , , , 

^,.(0) = ^^, CJpp{Q) = ^-, CTpqiO) = 0. (82) 

By using the complete positivity property of the dynamical semigroup $f, it was 
shown in Ref. ||2B| that the relation 

Tr{Mp)T.VjV,) = J2TriMp)V,^n^iMm) (83) 

j j 

represents the necessary and sufficient condition for p{t) = ^t{p) to be a pure state for 



all times t > 0. This equality is a generalization of the pure state condition |8^, ^, ^A 
to all Markovian master equations @). If p^(t) = p{t) for all t > 0, there exists a wave 
function ip &7i which satisfies the nonlinear Schrodinger type equation 

^^^ = H'm, (84) 



with the non-Hermitian Hamiltonian [^ 

H' = H+zY.< m.v.^m >^3-\< m^Y^yF^t) > ~lT.yIVr (85) 
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For environment operators Vj of the form (|^), the pure state condition (^) takes the 
following form ||2^, corresponding to equality in the relation (pO]): 

Dppaqq{t) + Dqqapp{t) - 2Dpg(rpg{t) = — (86) 

and the Hamiltonian (35) becomes 



*^—2 n //' , U\\2 , _ u\\ n //^ _/^N\2 



H' = H + \{ap{t)q - ag{t)f)) + -[A/i^ - Dpp{{q - ^^(t))^ + a,,(t)) - Z^,,((p - ^^(t)) 

+app(t)) + Dp,{{q - a,(t))(p - ap{t)) + (p - ap(t))(g - a,(t)) + 2ap,{t))]. (87) 

We will find the Gaussian states which remain pure during the evolution of the 



system for all times t. We start by considering the pure state condition (|86D and the 
generalized uncertainty relation ([68| ) which transforms into the following minimum 
uncertainty equality for pure states: 

(^,,it)app{t) - a^it) = -. (88) 

By eliminating app between the equalities ( P^j ) and (|5BD, like in Ref. |^, we obtain: 



I'^otI'^J n '' "^ n2 [I'^pgl'^J r) n - m ' 

^ \^ppJ-^qq ^pqJ 

Since the opening coefficients satisfy the inequality (^, we obtain from Eq. (|89|) the 
following relations which have to be fulfilled at any moment of time: 

DppD,q-Dl^ = ^, (90) 



h'^X 

Dppaqq{t) - Dpqapq{t) ^ = 0' (^l) 



h'^X 

apq{t){DppDqq - D'J - —Dpq = 0. (92) 
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From the relations ( ]88| ) and ( ]90| ) - ( p^ ) it follows that the pure states remain pure for 
all times only if their variances are constant in time and have the form: 

M^) = ^, M^) = ^, ^mW = ^. (93) 



If these relations are fulfilled, then the inequalities (|80D and (|8lD are both equivalent 
to (^, including also the corresponding equalities (^), ( FSl) and (pO|) . From Eq. 
(^) it follows that the variances remain constant and do not depend on time only if 
X(0) = X(oo), that means o"qg(0) = aqq{oo), crpp(O) = app{oo) ,apg{0) = apg{oo). Using 
the asymptotic values (^^ - (^61) of the variances and the relations (p3|) , we obtain 
the following expressions of the diffusion coefficients which assure that the initial pure 
states remain pure for any t : 

n =1^ D = ^^^^' D =-^ (94) 

Formulas (^^ are generalized Einstein relations and represent examples of quantum 
fluctuation-dissipation relations, connecting the diffusion with both Planck's constant 
and damping constant. With the coefficients (P^, the variances (123) - (EBI) become: 



""''^2^' """"^^kT' '''"^~2Q- ^ ' 



Therefore, the quantity a (see Eqs. (pOD , (|55D ) is equal to its minimum possible value 
^^/4, according to the generalized uncertainty relation (|79D . Then the corresponding 
state described by a Gaussian Wigner function is a pure quantum state, namely a 
correlated coherent state (squeezed coherent state) with the correlation coefficient 
([70|) r = —fi/u. Given cr^g, app and apg, there exists one and only one such a state 
minimizing the uncertainty a [Q. We remark that the minimization of the uncertainty 
a is equivalent, by virtue of the relations (^), ( ^5]) to the minimization of the entropy 
S. A particular case of our result (corresponding to A = /i and Dpg = 0) was obtained 
by Halliwell and Zoupas by using the quantum state diffusion method |^3|. We consider 
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here general coefficients /i and A and in this respect our expressions for the diffusion 
coefficients and variances generahze also the ones obtained by Dekker and Valsakumar 
\E^ and Dodonov and Man'ko B^, who used models where X = fi was chosen. If 



/i = 0, we get from ( P^ Dpq = 0. This case, which was considered in Ref. [^ and 

where we have obtained a density operator describing a pure state for any t, is also a 

particular case of the present results which give the most general Gaussian pure state 

which remains pure for any t. For /i = 0, the expressions (^) become 

h hmuj 

^m = ^^> ^PP = ^— ' ^pg = 0> (96) 

which are the values of variances (^) for the ground state of the harmonic oscillator 
and the correlation coefficient ( [70|) takes the value r = 0, corresponding to the case of 
usual coherent states. 

The Lindblad equation or its equivalent Fokker-Planck equation for the Wigner 
function with the diffusion coefficients (0) can be used only in the underdamped case, 
when oj > jji. Indeed, for the coefficients (0) the fundamental constraint (|^) implies 
that m'^{uj'^ — jj?)D'l > fi^X^ /A^ which is satisfied only if cj > /i. It can be shown that 
there exist diffusion coefficients which satisfy the condition (|90|) and make sense for 
uj < fx, but in this overdamped case we have always a > fi /A and the state of the 
oscillator cannot be pure for any diffusion coefficients [p7| . 

The fluctuation energy of the open harmonic oscillator is 

E{t) = ^^ppit) + -muj^a,,it) + fiapgit). (97) 



If the state remains pure in time, then the variances are given by ( pB]) and the fluctu- 
ation energy is also constant in time and is given by 

E = ji^Dpp + -mu^Dgg + fiDpg). (98) 



Minimizing this expression with the condition (|90|), we obtain just the diffusion coef- 



ficients (p^ ) and Emin = hVt/2. Therefore, the conservation of purity of state implies 
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that the fluctuation energy of the system has all the time the minimum possible value 
Emin- If the asymptotic state is a Gibbs state (H), we see from Eqs. (^ and (0) that 
the case when the diffusion coefficients satisfy the equality ( pOD corresponds to a zero 
temperature of the thermal bath and then the influence on the oscillator is minimal. 
In this limiting case ^ = and then Emm = huj/2, the correlation coefficient ([70|) van- 
ishes and therefore the correlated coherent state (squeezed coherent state) becomes 
the usual coherent state. 

If we choose the coefficients of the form (0), then the equation for the density 
operator can be represented in the form (^) with only one operator V, which up to a 
phase factor can be written in the form: 



with [V, ft] = 2n\. 

The correlated coherent states (^) with nonvanishing momentum average, can also 
be written in the form: 

1 1 1 2i i 

^^^^ " ^2^^~' ^""Pt-i^^^ " 7-^M)(a^ - ^,? + l^pA (100) 

and the most general form of Gaussian density matrices compatible with the general- 
ized uncertainty relation (|68D is given by Eq. (pj). For these matrices we can verify 
that Trp^ = h/2y/a and they correspond to the correlated coherent states (|100|) if 
o'qq, <Jpp and (Tpg in Eq. (^l]) satisfy the equality in Eq. (|68D . 

Consider the harmonic oscillator initially in a correlated coherent state (squeezed 
coherent state) of the form (|100|) , with the corresponding Wigner function (ffSl). For 



an environment described by the diffusion coefficients (p^), the solution for the Wigner 
function at time t is given by (pOf), where (Tq{t) and o"p(t) are given respectively by (|T6|) 
and (p!7D and the variances by (R^). Using Eq. (ffl|), we get for the density matrix the 
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following time evolution: 



,„,,, ,mQ.i , mQ,x + y ,,,9 



Ah' ''' h ' 2 "' ''' ^' h 

In the long time limit <Jq{t) = 0, crp{t) = and we have 

< x\p{oo)\y >= (^)^ exp{-^[f](x2 + y^) + z^(a;2 - y^)]}. (102) 



The corresponding Wigner function ( p2D has the form 

Woo{q,p) = ^exp[-— -(— / + mcj2^2 ^ 2/^^^)]. (103) 

7Tn Ml m 

We see that the time evolution of the initial correlated coherent state of the damped 

harmonic oscillator is given by a Gaussian density matrix with variances constant in 

time. According to known general results p9| , P0[ , the initial Gaussian density matrix 

remains Gaussian centered around the classical path. So, the correlated coherent state 

(squeezed coherent state) remains a correlated coherent state and cTqit) and crp{t) give 

the average time dependent location of the system along its trajectory in phase space. 

6 Decoherence and transition from quantum to clas- 
sical 

The environment induced quantum decoherence (loss of quantum coherence) and the 
transition from quantum to classical mechanics have recently been investigated in 
widely different contexts @, |, |, |, ^, |T|, P, ||, |19|, |0|, |2l|, ||, gSHI, H, |4|]. The 



central idea of the quantum decoherence approach is that the transition quantum - 
classical is a dynamical effect within quantum mechanics, keeping h ^ 0. The charac- 
teristic feature of the decoherence process is that an arbitrarily chosen generic initial 
quantum state will be affected on a characteristic decoherence time scale and only 
certain stable states (which turn out to be, in a sense, decay products of the other 
states) will survive during the time evolution of the system. 
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In quantum optics the quantum interference of Gaussian states in phase space 



was investigated by Schleich and coworkers p8| , p9|] and the decay of these quantum 
coherences was studied in Refs. [^, ^ (for a recent review on quantum interference 
in phase space see the work of Buzek and Knight 



The predictabihty sieve was recently proposed ^0\ and implemented for a harmonic 



oscillator with the resulting evolution of the reduced density matrix [T^, 19|. For a 



weakly damped harmonic oscillator, pure states selected by the predictability sieve, 



called preferred states, turn out to be the coherent states |^, ^. Decoherence is 
caused by the loss of phase coherence between the preferred quantum states in the 
Hilbert space of the system due to the interaction with the environment 0, ^. Pre- 
ferred states are singled out by their stability (measured, for example, by the rate of 
predictability loss - the rate of entropy increase) under the joint influence of the envi- 



ronment and the self-Hamiltonian [21]. Thus, the strength and nature of the coupling 
with the environment play a crucial role in selecting preferred states, which - given 
the distance-dependent nature of typical interactions - explains the special function 
of the position observable [Q, ^. Coupling with the environment also sets the deco- 
herence time scale - the time on which quantum interference between preferred states 
disappears p|, §, |l^, ^, ^, ^ . Classicality is then an emergent property of an open 
quantum system. 

The interaction with the environment also induces fluctuations in the evolution of 
the system. There is therefore a certain degree of tension between the demands of 
decoherence and approximate classical predictability: decoherence requires interaction 
with an environment, which inevitably produces fluctuations, but classical predictabil- 
ity requires that these fluctuations be small. The von Neumann entropy is an ideal 
parameter with which to characterize the decoherence, i. e. the rapid decay of off- 
diagonal coherences. In [^ an information-theoretic measure of the size of these 
fluctuations was proposed - the Shannon information of the Qh distribution of the 
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density matrix p for the distinguished system (see Sec. 4). This measure of uncertainty 
is in fact bounded from below by the von Neumann entropy of p. This suggests that 
the von Neumann entropy is the key to understanding the connection between deco- 
herence and fluctuations: it limits the amount of decoherence from above but bounds 
the size of the fluctuations from below. Large entropy, which is a signal of destruction 
of interference, therefore permits good decoherence, but leads to large fluctuations; on 
the other hand, small entropy allows small fluctuations, but the amount of decoherence 
is also small. Environmentally induced fluctuations are inescapable if one is to have 
decoherence. Although there is some tension, there is a broad compromise regime in 
which decoherence and classical predictability can each hold extremelly well ||4CI|| . 

As we have already mentioned, besides the von Neumann entropy S (^), there is 
another quantity which can measure the degree of mixing or purity of quantum states. 
It is the linear entropy Si (|63|). For pure states Si = Q and for a statistical mixture 
Si> 0. The increasing of the linear entropy Si (as well as of the von Neumann entropy 
5*) due to the interaction with the environment is associated with the decoherence 
phenomenon, given by the diffusion process ||T^ 120]. Due to the decoherence, the 



macroscopic systems obey essentially classical equations of motion, despite the quan- 
tum mechanical nature of the underlying microscopic dynamics. Dissipation increases 
the entropy and the pure states are converted into mixed states. The rate of entropy 
production is given by 

Slip) = -2Tr(p^) = -2Tr(pL(p)), (104) 

where L is the evolution operator. According to Zurek's theory, the maximally predic- 
tive states are the pure states which minimize the entropy production in time. These 
states remain least affected by the openness of the system and form the preferred set of 
states in the Hilbert space of the system, known as the pointer basis. Their evolution 
is predictible with the principle of least possible entropy production. 
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Using Eq. (|]), in our model the rate of entropy production ( |104| ) is given by: 

4 
Slit) = ■^[DppTi{p^'^ - pqpq) 

+^ggTr(/3V - PppP) - DpgTi{p^{qp + pq) - 2pqpp) - —Tt{p'^)]. (105) 

When the state remains approximately pure (p^ ~ p), we obtain: 

Slit) = -2 (Dpp(T,,(t) + Dggappit) - 2Dpgapg{t) - — -) > 0, (106) 



according to (|80D . We see that Siit) = when the condition (|8q) of purity for any time 
t is fulfilled. The entropy production Si is also equal to at t = if the initial state is 
a pure state. We have shown in the preceding section that the only initial states which 
remain pure for any t are the correlated coherent states (squeezed coherent states) and 
therefore, we can state that in the Lindblad theory for the open quantum harmonic 
oscillator the correlated coherent states, which are generalized coherent states, are the 
maximally predictive states. Our result generalizes the previous results which assert 
that for many models of quantum Brownian motion in the high temperature limit the 
usual coherent states correspond to minimal entropy production and, therefore, they 
are the maximally predictive states. In our model the coherent states can be obtained 
as a particular case of the correlated coherent states by taking p = 0, so that the 
correlation coefficient r = (see Eq. (|70|) ). 

Paz, Habib and Zurek |l^, Q considered the harmonic oscillator undergoing quan- 
tum Brownian motion in the Caldeira-Leggett model and concluded that the minimiz- 
ing states which are the initial states generating the least amount of von Neumann or 
linear entropy and, therefore, the most predictible or stable ones under evolution in the 
presence of an environment are the ordinary coherent states. Using an information- 
theoretic measure of uncertainty for quantum systems, Anderson and Halliwell showed 



in ||39[ that the minimizing states are more general Gaussian states. Anastopoulos and 



Halliwell ^^ offered an alternative characterization of these states by noting that these 
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states minimize the generalized uncertainty relation. According to this assertion, we 
can say that in the Lindblad model the correlated coherent states (squeezed coherent 
states) are the most stable states which minimize the generalized uncertainty relation 
(|55|). Our result confirms that one of |^, but the model used in |^ is different, 



namely the open quantum system is a particle moving in a harmonic oscillator poten- 
tial and linearly coupled to an environment consisting of a bath of harmonic oscillators 
in a thermal state. At the same time we remind that the Caldeira-Leggett model con- 
sidered in [0, ^ violates the positivity of the density operator at short time scales 
||53| , Q, whereas in the Lindblad model considered here the property of positivity is 
automatically fulfilled. 

The rate of predictability loss, measured by the rate of linear entropy increase, for 
a damped harmonic oscillator is also calculated in the framework of Lindblad theory 
in Ref. ^5\. The initial states which minimize the predictability loss are identified as 



quasi-free states with a symmetry dictated by the environment diffusion coefficients. 
For an isotropic diffusion in phase space, the coherent states or mixtures of coherent 
states are selected as the most stable ones. 

In order to generalize the results of Zurek and collaborators, the entropy production 
was considered by Gallis [Q within the Lindblad theory of open quantum systems, 
treating environment effects perturbatively. Gallis considered the particular case with 
Dpq = and found out that the squeezed states emerge as the most stable states for 
intermediate times compared to the dynamical time scales. The amount of squeezing 
decreases with time, so that the coherent states are most stable for large time scales. 
For Dpq 7^ our results generalize the ones of Gallis and establish that the corre- 
lated coherent states are the most stable under the evolution in the presence of an 
environment. 
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7 Summary and concluding remarks 

Recently there is a revival of interest in quantum Brownian motion as a paradigm of 
quantum open systems. The possibility of preparing systems in macroscopic quantum 
states led to the problems of dissipation in tunneling and of loss of quantum coher- 
ence (decoherence). These problems are intimately related to the issue of quantum 
to classical transition and all of them point the necessity of a better understanding 
of open quantum systems. The Lindblad theory provides a selfconsistent treatment 
of damping as a general extension of quantum mechanics to open systems and gives 
the possibility to extend the model of quantum Brownian motion. In the present pa- 
per we have studied the one-dimensional harmonic oscillator with dissipation within 
the framework of this theory. From the master equation of the damped quantum 
oscillator we have derived the corresponding Fokker-Planck equation in the Wigner 
W representation. The obtained equation describes an Ornstein-Uhlenbeck process 
and the Wigner function is a two-dimensional Gaussian. We have then obtained the 
explicit form of the density operator from the master and Fokker-Planck equations. 
The density operator in a Gaussian form is a function of the position and momentum 
operators in addition to several time dependent factors. In the long time limit the 
density operator approaches a thermal state. Then the density operator has been used 
to calculate the von Neumann entropy and the effective temperature. The temporal 
behaviour of these quantities shows how they approach their equilibrium values. Fol- 
lowing Anderson, Anastopoulos and Halliwell |]39| , ^ , we have put the von Neumann 
entropy in association with uncertainty and linear entropy. The Wehrl entropy gener- 
ates an informatical-theoretical measure of the size of the intrinsic state fluctuations. 



As coherent states are most robust in dissipative environment |T^, [2^, ^, ^, this 
suggests the utility of the Wehrl entropy in characterizing the decoherence process in 
quantum mechanics. We have also shown that the only states which stay pure during 
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the evolution in time of the system are the correlated coherent states (squeezed co- 
herent states) under the condition of a special choice of the environment coefficients. 
These states are also connected with the decoherence phenomenon and they are the 
most stable under the evolution in the presence of the environment. The obtained 
results in the framework of the Lindblad theory can be used for the description in 
more details of the connection between uncertainty, decoherence and correlations of 
open quantum systems with their environment. 
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